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Abstract 

Using thermodynamic arguments we find that the probability that there are no 
eigenvalues in the interval (— s, oo) in the double scaling limit of Hermitean matrix 
models is O (exp(— s 27+2 )) as s — > +oo. Here 7 = m — 1/2, m = 1, 2, • • • determine 
the m th multi-critical point of the level density: cr(x) oc b [1 — (x/6) 2 ] 7 , x 6 (—6,6), 

6 2 oc JV. Furthermore, the size of the transition zone where the eigenvalue density 

7-1 

becomes vanishingly small at the tail of the spectrum is ~ jV 2 ^ +1 ) in agreement 
with earlier work based upon the string equation. 

To appear as a Letter in J. Phys. A: Math. Gen. 1995. 
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A basic quantity in Hermitean matrix models is the probability, £^(0; J), that a set J 
contains no eigenvalues. For N x N Hermitean matrix models with unitary symmetry we 
have the well-known expression 

with 

dfi(x) = Yl \x a -x b \ 2 Yl dx a , 

\<a<b<N l<a<N 

J the complement of J and V(x) is the "confining" potential [M. As indicated in (1), 
minus the logarithm of this probability has the physical interpretation, in terms of Dyson's 
Coulomb fluid || 0, [12|, || , as the change in free energy 



AF = F[J] — F[J U J]; (2) 

that is, the free energy of the iV charges confined to region J minus the free energy of iV 
charges in the natural support J U J of w(x) := e~ v ^ x \ 

In this paper we shall mainly consider the case J = (— s, oo), s > 0, and write E 2 (s) for 
-^2(0; (— s, oo)). Below we shall use the continuum approximation of Dyson which treats 
the iV eigenvalues in the large iV limit as a continuous fluid described by a continuous 
charge density a with the free energy expressed in terms of a 0. This approximation has 
been previously applied to the Unitary Laguerre Ensemble (ULE) where w(x) = x a e~ x , 



2 



X 



G (0, oo) and a > — 1 [[]].[] Here we examine matrix models witl 



with g2 — 1. In principle we should not have to make the continuum approximation since 
it is known that ^(O; J) is expressible in terms of solutions to a completely integrable 



system of partial differential equations |19| . However, the analysis of these equations for 
V of the above form is quite difficult. (Of course, the Gaussian case is not included in this 
remark.) It is hoped that the approximate expressions derived here, which we believe are 
asymptotic as s — > oo, will aid in the analyses of these equations. 

To begin, consider the Gaussian Unitary Ensemble (GUE) with #2fc+2 = for k > 1. 
For the scaled GUE with J = (— t, t) it is known that E 2 (0; {—t, t)) is a r-function of a par- 



ticular fifth Painleve transcendent [ I5| . Starting with this representation an asymptotic 
expansion for £^(0; (—£,£)) as t — * oo can be derived, though the first such asymptotic 
expansion was achieved by Dyson using methods of inverse scattering [11]. (Actually, 



1 It is known from the theory of liquids (by an application of the Boltzmann principle) that the 
probability, Pd(R), of finding a bubble of radius R in the bulk of a fluid (in d-dimensions) at equilibrium 
with temperature 1/(3 is 

P d (R) - e -PEvR d -f)E av R d - 1 ^ R > coherence length, 

where Ey is the energy/volume for creating a bubble and Egy is the surface energy. If we specialise this 
formula to d = 1 then 

Pi(R) _ e - constantfl , 

in contradiction with the known result This is due to the fact the Coulomb fluid has long ranged 
interactions. 

2 The reason for our choice of notation for the coefficients of V will become clear below. 
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there is still an undetermined constant from either the inverse scattering analysis or the 
Painleve analysis, see, e.g. 0.) The leading term, — In E 2 (0] (— t, t j) ~ ^7r 2 t 2 , was first 
obtained from the fluid approximation |9|]. Indeed, the t 2 term of the asymptotic expan- 
sion can be given a simple physical interpretation: it is proportional to the square of the 
number of eigenvalues excluded in the (scaled) interval (— t, t) since in the bulk scaling 
limit of the GUE the eigenvalue density is a constant ~ xML_ This suggests that a nat- 
ural variable is one which gives uniform density in the excluded interval. We can always 
achieve this by a simple change of variables since the problem is one-dimensional. By 
introducing a new variable £ and a corresponding via the relation 

f?(Odf := 1 • <% = a(x)d Xj (4) 
the density in the £ "scale" is made unity. Therefore, —ha.E 2 (0;J) is asymptotic to 

2 

, J = (xx,x 2 ). 

We conclude from the above arguments that for a large interval, 

— \nE 2 (0; J) ~ N 2 (l), (5) 

where 

N(l) = number of eigenvalues excluded in an interval of length /. (6) 

We mention that a screening theory of the continuum Coulomb fluid gives a physical 
justification of these arguments [0] though we know of no proof of the general validity of 
this relationship. 

To further test the validity of (Ejl) and ([5] ) we consider the edge scaling limit of the 



/ d£ 


2 


' r%2 




/ a(x) dx 









GUE where exact results are known |T7 1 . Accordingly, we simply compute the number of 
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eigenvalues excluded from an interval of length I (= b — a) from the soft edge, b = V2N, 

3.1 3 / 2 3/2 



iV(Z) = / dx-Vb 2 -x 2 oc v 7 ^ / dxVt^x oc [2 1/2 iV 1/6 Z 

Ja IT Ja L 

Observe that N 2 (l) ~ s 3 not only supplies the correct exponent in the scaled variable s (= 
2 1 / 2 iV 1//6 /) in — lni? 2 (s), we also have the correct density at the soft edge: a^{-\/2N) = 
2 i/2jvV6 j w hich agrees with known exact results |16l O. This result predicts the shrinking 



of the size of the transition zone (~ N- l ' % ) as N — > oo — a reasonable behaviour from the 
Coulomb fluid point of view since the GUE potential x 2 is strongly confining. The same 
approximation has been applied to the origin scaling limit of the ULE H and the result 



agrees with the first term of the exact asymptotic expansion [IS . 

These two confirmations of the validity of and (0) give us confidence to apply the 
method to the matrix models with V given by (H). (These are the cases of interest in the 



matrix models of 2D quantum gravity [^|, The charge density a satisfies an integral 

equation||, [TJJ derived from the following minimum principle: 

min a F[a], 

F[a] = J dxV(x)a(x) — J dx J dya(x) In \x — y\cr(y) (7) 
subject to the constraint Jj dxa(x) = N, which is 

rb 

V(x) — 2 / dyla \x — y\cr(y) = constant = chemical potential, x E (—6, b). (8) 

J ~ b 

Since V is even so is a, and making use of this symmetry (P) becomes 

V(x) — 2 / dy In \x 2 — y 2 \cr(y) = constant. (9) 
Jo 

With the change of variables x 2 = u, y 2 = v and r(u) = a(y/u)/(2y^u), @ becomes 

V(\/u) — 2 / dv r(v) In \u — v\ = constant, uG(0,6 2 ). (10) 
Jo 
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This is converted into a singular integral equation by differentiating with respect to u: 

dv ^ 2Pr 2 <fo^i = o, « e( o,i"). (11! 



du Jo u — v 

Here b, which determines the upper and lower band edges, is fixed by the normalization 
condition J b b a(x)dx = N . 

Following m the solution isQ 



1 Jb 2 -u^f b2 dv I v dV(yfc) 2N 



2tt z V u Jo v — u V tr — u em 



'b 2 — u / , 3 u \ 

5:^2^1 , (12) 



« H ' '2' 62 



where the integral can be found in [13] and 



Returning to a, it can be shown that 



X\ 2 



X N 2 



(13) 



where n p (z) is a polynomial of degree p in 2 with coefficients depending on the linear 
combinations of the coupling constants g^. The edge parameter b is determined from the 
normalization condition and reads b 2 = CN where 

c 



dtVT^t 2 U p {t 2 ) 
is independent of N. 



3 constant so i ves ^ e homogeneous part of (ITT|) . However, based on the variational principle, including 

^/u(b 2 — u) 1 — v 

this solution would increase the free energy. 2-Pi( — k, 1, 3/2, z) = X)n=o (3/2)" z ™ ^ s a polynomial of degree 
k in 2. 



Taking the special case p = 1 (now g± = g), we have 
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By tuning g to g c , such that — g c = 1 + g c /2, we have 



a(x) = constant b 



x 



3/2 



producing a qualitative deviation in the density at the edges (±6) of the spectrum from 
the Wigner's semi-circle distribution || Q. A calculation now gives 



N{1) CC f'dx b(l - x/bf'\l + x/bfl* » A f\ h _ x) 3/2 „ ( _L_ 



5/2 



:s 5 / 2 , (14) 



and thus — In E 2 (l) ~ s 5 . Observe that due to this tuning the length of the transition 
zone (~ is now an increasing function of N. It is clear that the tuning procedure 

can be generalized to p > 10. By simultaneously adjusting the coupling constants g 4 , g 6 



etc., to their respective critical values we can have 



a(x) = constant b 



1 - 



x 



21 T 



(15) 



where 7 = P + Computing iV(Z) we find 
N(l) oc J^dxbU 



I 



-Y (i + -V oc ( 



7+1 



7+1 



(16) 



Therefore 



log E 2 (s) « -s 



27+2 



S — >■ CO 



(17) 



The non-perturbative soft edge density is determined as 



1-7 



a N (VN) » jV^a+^T, ^ oo. 



In the quantum gravity literature 7 = m — 1/2, m = 1, 2 • 



The corresponding size of the transition zone is ~ iW, where 

7-1 
^"2(7+1) 

a result previously obtained from the string equation ]5], |8|]. Note that our x variable is 
related to Bowick and Brezin's ||] A as x — y/NX. Supplying the appropriate y/N factor 



we obtain from (IS) Bowick and Brezin's result N 2m+i. 
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